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Abstract. We study dispersive properties for the wave equation in the Schwarzschild 
space-time. The first result we obtain is a local energy estimate. This is then used, 
following the spirit of |29], to establish global-in-time Strichartz estimates. A consid- 
erable part of the paper is devoted to a precise analysis of solutions near the trapping 
region, namely the photon sphere. 



1. Introduction 



The aim of this article is to contribute to the understanding of the global-in-time 
dispersive properties of solutions to wave equations on Schwarzschild black hole back- 
grounds. Precisely, we consider two robust ways to measure dispersion, namely the local 
energy estimates and the Strichartz estimates. 

Let us begin with the local energy estimates. For solutions to the constant coefficient 
wave equation in 3 -|- 1 dimensions, 

Du—0, u{0) = Uq, Ut(0) = Ml, 

we have the original estimates of Morawetz |33Q 

(1.1) /7 l-\y7u\^{t,x)dtdx<\\Vuo\\h + \\u^\\l, 
Jo fI 

where y/ denotes the angular derivative. To prove this one multiplies the wave equation 
by the multiplier {dr + and integrates by parts. Within dyadic spatial regions one 
can also control u, dtu and drU. Precisely, we have the local energy estimates 

(1.2) i?"^||Vu||i2(jjxB(o,fl)) + R~^\W\\l^(R-kB(o,H)) ^ I|Vuo||l2 + 
See for instance [20], [22], [iQ], [M], [12]. 

One can also consider the inhomogeneous problem, 

(1.3) = u(0)=wo, Mt(0)=Mi, 

The authors were supported in part by the NSF grants DMS0354539 and DMS0301122 . 

^There is another estimate commonly referred to as a Morawetz estimate. This corresponds to using 
the multiplier (t^ -I- r^)9t -I- 2trdr. We will reserve the term Morawetz estimate for jl.l) and shall call 
the latter estimate the Morawetz conformal estimate. 
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In view of (1.2| we define the local energy space LEm for the solution u by 



(1-4) \\u\\le,, ^sup\2 3||Vu|U2(^^) + 2 '^\\u\\l2(^a,) 

where 

= M X {2^ < \x\ < 2^+1}.. 
For the inhomogeneous term / we introduce a dual type norm 

\\f\\LEl,^J2^'^\\f\\LHA,y 

Then we have: 



Theorem 1.1. The solution u to (1.3 I satisfies the following estimate: 
(1-5) Mle,, < llVuolU^ + II^iIIl^ + II/IIlb-, 

One may ask whether similar bounds also hold for perturbations of the Minkowski 
space-time. Indeed, in the case of small long range perturbations the same bounds as 
above were established very recently by two of the authors, see |301 Proposition 2.2] 
or [28^ (2.23)] (with no obstacle, = 0). See also [1], for related local energy 
estimates for small perturbations of the d'Alembertian. For large perturbations one faces 
additional difficulties, due on one hand to trapping for large frequencies and on the other 
hand to eigenvalues and resonances for low frequencies. The Schwarzschild space-time, 
considered in the present paper, is a very interesting example of a large perturbation of 
the Minkowski space-time, where trapping causes significant difficulties. 

The Schwarzschild space-time is a spherically symmetric solution to Einstein's 
equations with an additional Killing vector field K, which models the exterior of a massive 
spherically symmetric body. Factoring out the component it can be represented via 
the Penrose diagram: 

The radius r of the spheres is intrinsically determined and is a smooth function on 
A4 which has a single critical point at the center. The regions / and /' represent the 
exterior of the black hole, respectively its symmetric twin, and are characterized by the 
relation r > 2M. We can represent / as 

/ = M X (2Af , oo) X §2 

with a metric whose line element is 



2 



(1.6) ds'^-{l-'^yt^+(l-'^) dr' + r'du; 

where duj"^ is the measure on the sphere S^. The Killing vector field K is given by 
K = dt, which is time- like within /. The differential dt is intrinsic, but the function t is 
only defined up to translations on /. 

The regions // and //' represent the black hole, respectively its symmetric twin, the 



white hole, and are characterized by the relation r < 2M. The same metric as in (1.6 1 
can be used. The Killing vector field K is still given hy K = dt, which is now space- like. 
Light rays can enter the black hole but not leave it. By symmetry light rays can leave 
the white hole but not enter it. 
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Figure 1. The Penrose diagram for the Kruskal extension of the 
Schwarzschild solution 

The surface r = 2M is called the event horizon. While the singularity at r = is 
a true metric singularity, we note that the apparent singularity at r = 2M is merely a 
coordinate singularity. Indeed, denote 

r* =r + 2M log(r - 2M) - 3M - 2M log M, 

so that 

dr* = (^1 j dr, r*(3M) = 

and set V = t + r* . Then in the {r,v,oj) coordinates the metric in region / is expressed 
in the form 

ds^ = - (l - '^)dv^ + 2dvdr + r^dw^ 

which extends analytically into the black hole region 1 + II. In particular, given a choice 
of the function t in region /, this uniquely determines the function t in the region II via 
the same change of coordinates. 

In a symmetric fashion we set w = t — r*. Then in the (r, w, uj) coordinates the metric 
is expressed in the form 

ds^ = - (l - ^)'^'*"^ - 2dwdr + r'^dJ^, 
which extends analytically into the white hole region 7+ 

One can also introduce global nonsingular coordinates by rewriting the metric in the 
Kruskal-Szekeres coordinate system, 

V =e*", w = — e . 
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However, this is of less interest for our purposes here. Further information on the 
Schwarzschild space can be found in a number of excellent texts. We refer the inter- 
ested reader to, e.g., [18], [31], and |51] . 

As far as the results in this paper are concerned, for large r the Schwarzschild space- 
time can be viewed as a small perturbation of the Minkowski space-time. The difficulties 
in our analysis are caused by the dynamics for small r, where trapping occurs. The 
presence of trapped rays, i.e. rays which do not escape either to infinity or to the 
singularity r = 0, are known to be a significant obstacle to proving local energy, dispersive, 
and Strichartz estimates and, in some case, are known to necessitate a loss of regularity. 
See, e.g., [10] and [57] . 

There are two places where trapping occurs on the Schwarzschild manifold. The first is 
the surface r = 3M which is called the photon sphere. Null geodesies which are initially 
tangent to the photon sphere will remain on the surface for all times. Microlocally the 
energy is preserved near such periodic orbits. However what allows for local energy 
estimates near the photon sphere is the fact that these periodic orbits are hyperbolic. 
The second is at the event horizon r — 2M, where the trapped geodesies are the vertical 
ones in the {r,v,uj) coordinates. However, this second family of trapped rays turns out 
to cause no difficulty in the decay estimates since in the high frequency limit the energy 
decays exponentially along it as u ^ oo. This is due to the fact that the frequency decays 
exponentially along the Hamilton flow, and in the physics literature it is well-known as 
the red shift effect. 

To describe the decay properties of solutions to the wave equation in the Schwarzschild 
space, it is convenient to use coordinates which make good use of the Killing vector fleld 
and are nonsingular along the event horizon. The (r, v, cu) coordinates would satisfy these 
requirements. However the level sets of v are null surfaces, which would cause some minor 
difficulties. This is why in / -I- // we introduce the function v defined by 

V — V — ii(r) 

where /j, is a smooth function of r. In the {v,r,uj) coordinates the metric has the form 

ds' = -(l - ?M)d^)2 ^2 - (l - ^)m'(^)) dvdr 

+ (2^'(r) - (l - ^) dr' + r^dLo\ 

On the function /i we impose the following two conditions: 

(i) /i(r) > r* for r > 2M, with equality for r > 5M /2. 

(ii) The surfaces v — const are space-like, i.e. 

//'(r)>0, 2-{l-'^)^^\r)>0. 

The first condition (i) insures that the {r,v,uj) coordinates coincide with the {r,t,uj) 
coordinates in r > 5M/2. This is convenient but not required for any of our results. 
What is important is that in these coordinates the metric is asymptotically flat as r ^ oo. 

I„ .he p™t 0, .he S.riohart. i. . al» .ha. = (l - ^1 " nea. 
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r = 3M, which in other words says that we can work in the (r, t) coordinates near the 
photon sphere. However, this may be merely an artifact of our method. 

We introduce a symmetric function vi in /' + II, as well as the the functions w and 
fill in / + //', respectively /' + //'. Given a parameter < rp < 2M we partition the 
Schwarzschild space into seven regions 

M = MrU Ml^ M'nU M'lU Mt^ Mc^ Mb 
as in Figure [2] The right/left top/bottom regions are 

r = ro<2M 



= 3M V 
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r-Q < 2M 

Figure 2. The Schwarzschild space partition represented on the Pen- 
rose diagram 



Mb. = {v>0, r > ro} C / + //, Ml = {wi > 0, r > ro} C /' + //, 

M'ji = {w<0, r>ro}cI + II', M'^ = {wi<0, r> ro} C /' + //', 
the top and bottom regions are 

Mt ^ {r < ro} nil, Mb ^ {r < ro} n //', 
and the central region Mc '^s the remainder of M. Moreover, define 

= Mr n{v = 0} 

E+^MRn{r^ro}. 
and similarly for the other regions. 

In what follows we consider the Cauchy problem 

(1-7) Ogc^^f, 0|so=0o, K<Pi^,^(j>i 

where for convenience we choose the initial surface Eq to be the horizontal surface of 
symmetry 

Eo = {t-o}n(/ + /') 
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and K is smooth, everywhere timehke and equals K on Sq outside Aic- Observe that 
we cannot use K on all of Sq since it is degenerate at the center (i.e. on the bifurcate 
sphere). 



The equation (1.7 1 can be solved as follows: 

(i) Solve the equation in Mc with Cauchy data on Sq. Since A^c is compact and has 
forward and backward space-like boundaries, this is a purely local problem. 

(ii) Solve the equation in Aifi with Cauchy data on T,~^. The forward boundary of 

M 

R is S^, which is space-like. This is the most interesting part, where we are interested 
in the decay properties as v ^ oo. In a similar manner solve the equation in AIl, M.'j^ 
and TW'^. 

(iii) Solve the equation in Mt with initial data on the space-like surface — ^ 
To} n II. Here one can track the solution up to the singularity and encounter a mix of 
local and global features. This part of the analysis in not pursued in the present article. 

A significant role in our analysis is played by the Killing vector field K, which in the 
(r, 5) coordinates equals This is time-like outside the black hole but space-like inside 
it. Furthermore, it is degenerate at the center. Using the Killing vector field outside 
the black hole we obtain a conserved energy i?o[</>] for solutions (f> to the homogeneous 
equation Dg^ = 0. On surfaces t — const in the (r, t) coordinates the energy £'o[0](i) has 
the form 



(1.8) Eo[<j^]= I '(9*0)2+ (l-—)(a.</.)2 + 1)^01^ 



2 



drduj. 



Since the vector field K is degenerate at the center, so is the corresponding energy 
at r = 2M. Hence it would be natural to replace it with a nondegenerate energy, which 
on the initial surface Eq can be expressed as 



(1.9) i?H(So) = / / [(l - ^) ' (dtcbf + (l - ^) ' {dr^r + \f<P\ 



drduj. 



Unfortunately this is no longer conserved, and this is one of the difficulties which we face 
in our analysis. We remark that a related form of a nondegenerate energy expression was 
introduced in |T1] and proved to be bounded in the exterior region on surfaces t = const. 

Part of the novelty of our approach is to prove bounds not only in the exterior region, 
but also inside the event horizon. This is natural if one considers the fact that the 
singularity at r = 2M is merely a removable coordinate singularity. In order to do this, 
it is no longer suitable to measure the evolution of the energy on the surfaces t = const 
(see below). Thus the above energy i?[<^](So) is relegated to a secondary role here and is 
used only to measure the size of the initial data. 

A priori the energy E[(t)]{t) of only determines its Cauchy data at time t modulo 
constants. However, in what follows we implicitly assume that decays at oo, in which 
case (j) can be also estimated via a Hardy-type inequality. 



(1.10) J [1 j r^drduj ^ J j idr<P) r^drduo. 

This is proved in a standard manner; the details are left to the reader. 
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We shall now further describe our main estimates in the region AAr. the local energy 
decay, the WKB analysis which yields a local energy decay with only a logarithmic loss, 
and finally the Strichartz estimates. 

For the initial energy on S]^ we use 

E[cf\{^~^)^ f {\dr<l>\^ + \d,cl>\^ + \y7<l>\^)r^drdiu. 

For the final energy on we set 



We also track the energy on the space-like slices v — const, 

J MRr]{v=vo) 

Thus E[cj>]{ll-^) = E[4>m. 

For the local energy estimates one may first consider a direct analogue of the Minkowski 
bound (1.5 I. Unfortunately such a bound is hopeless due to the trapping which occurs 
at r = 3M. Instead, for our first result we define a weaker preliminary local energy space 
LEa with norm 
(1.11) 



\LEa 




+ + — r^drdvduj. 

r y r I 



Compared to the LEm norm we note the power loss in the angular and v derivatives 
at r = 3Af. The LEq norm is also weaker than LEm as r — > 00, but this is merely for 
convenience. 

At the same time we would like to also consider the inhomogeneous problem Dgcf) = f . 
To measure the inhomogeneous term /, we introduce the norm LEq, which is stronger 
than LE\.f 

(1.12) ll/lli = / (l-^-^Y r^f r'drdvdoj. 



Mi 



r 



Again the important difference is at r = 3Af. Our first local energy estimate is the 
following: 

Theorem 1.2. Let (j) solve the inhomogeneous wave equation Og(f) — f on the Schwarzs- 
child manifold. Then we have 

(1.13) Emi:+)+ sup E[4>m + \miE,<E[mB) + \\f\\iE5- 

v>0 

Here we made no effort to optimize the weights at r = 3Af and r — 00. This is 
done later in the paper. On the other hand the above estimate follows from a relatively 
simple application of the classical positive commutator method. The advantage of having 
even such a weaker estimate is that it is sufficient in order to allow localization near the 
interesting regions r = 3M and r = 00, which can then be studied in greater detail using 
specific tools. 
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The first related results regarding the solution of the wave equation on Schwarzschild 
backgrounds were obtained in [50j and [24| which proved uniform boundedness in region 
/ (including the event horizon). The first pointwise decay result (without, however, a 
rate of decay) was obtained in jl9]. Heuristics from j36] suggest that solutions to the 
wave equation in the Schwarzschild case should locally decay like v~^. For spherically 
symmetric data a decay rate was obtained in |l6j, and under the additional as- 

sumption of the initial data vanishing near the event horizon, the decay rate was 
proved in |23]. In general the best known decay rate, proved in [H], is (see also [3)- 
We also refer the reader to p^, where optimal pointwise decay rates for each spherical 
harmonic are established for a closely related problem. 



Estimates related to (1.131 were first proved in |25] for radially symmetric Schrodinger 
equations on Schwarzschild backgrounds. In [H |3] S], those estimates are extended to 
allow for general data for the wave equation. The same authors, in [316], have provided 
studies that give improved estimates near the photon sphere r — 3M. 

Moreover, we note that variants of these bounds have played an important role in 
the works [7] and [Tl] which prove analogues of the Morawetz conformal estimates on 
Schwarzschild backgrounds. This allows one to deduce a uniform decay rate for the local 
energy away from the event horizon, though there is necessarily a loss of regularity due to 
the trapping that occurs at the photon sphere. Instead in this paper we restrict ourselves 
to time translation invariant estimates, and we aim to clarify/streamline these as much 
as possible. 

All of the above articles use the conserved (degenerate) energy i?o[</>] on time slices, 
obtained using the Killing vector field dt ■ As such, their estimates are degenerate near 
the event horizon. Further progress was made in jlj, where an additional vector field 
was introduced near the event horizon, in connection to the red shift effect. This enabled 
them to obtain bounds in the exterior region involving a nondegenerate form of the energy 



related to (1.9 1 



The approach of |25], p], [7] and |T3] is to write the equation using the Regge- Wheeler 
tortoise coordinate and to expand in spherical harmonics. For the equation corresponding 
to each spherical harmonic, one uses a multiplier which changes sign at the critical point 
of the effective potential. 

Here we work in the coordinates (r, v, cu), though this is not of particular significance. 



and we do not expand into spherical harmonics. We prove (1.13 1 using a positive commu- 
tator argument which requires a single differential multiplier. We hope that this makes 
the methods more robust for other potential applications. 

During final preparations of this article, localized energy estimates proved without 
using the spherical harmonic decomposition also appeared in [15]. The methods contained 
therein are somewhat different from ours. 



Compared to the stronger norms LEm, LE^j the weights in (1.131 have a polynomial 
singularity at r = 3M, which corresponds to the family of trapped geodesies on the 
photon sphere. As a consequence of the results we prove later, see Theorem |3.2| the 
latter fact can be remedied to produce a stronger estimate. 

Theorem 1.3. Let (f> solve the inhomogeneous wave equation ^n<t> ^ f on the Schwarzs- 



child manifold. Then (1.131 still holds if the coefficient (1 — 3A//r) in the LEq and the 
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^ 3M 






) 




r 





LEq norms is replaced by 



Now we have only a logarithmic singularity at r = 3M. The result above is only stated 
in this form for the reader's convenience. The full result in Theorem |3.2| is stronger but 
also more complicated to state since it provides a more precise microlocal local energy 
estimate. 

The logarithmic loss is not surprising, since it is characteristic of geometries with 
trapped hyperbolic orbits (see for instance |9], [12], |3l]). Indeed, a similar estimate 
in the semiclassical setting is obtained in |13) using entirely different techniques. Note, 
however, that the aforementioned estimate only involves logarithmic loss of the frequency; 
our result is stronger since it also implies bounds for |j(ln |r*|)~^u||i2, which are necessary 
in order to prove Strichartz estimates. 

There are two regions on which the analysis is distinct. The metric is asymptotically 
flat, and thus, near infinity, one can retrieve the classical Morawetz type estimate. On 
the other hand, around the photon sphere r = 3M we take an expansion into spherical 
harmonics as well as a time Fourier transform. Then it remains to study an ordinary 
differential equation which is essentially similar to 

{dl-X'{x' + e))u^f, |x|<l. 

For this we use a rough WKB approximation in the hyperbolic region combined with 
energy estimates in the elliptic region. Airy type dynamics occur near the zeroes of the 
potential. 

Even though it is weaker, the initial bound in Theorem |1.2| plays a key role in the 
analysis. Precisely, it allows us to glue together the estimates in the two regions described 
above. 

We next consider the Strichartz estimates. For solutions to the constant coefficient 
wave equation on M x M'^, the well-known Strichartz estimates state that 

(1.14) \\\D,\-''^vu\\^.^^.^ < wvumiL^ + iii^xrvii^.i^,^ 

Here the exponents {pi,pi,qi) are subject to the scaling relation 

1 3 3 

and the dispersion relation 

1 1 1 

1.16 - + -<-, 2<p<oo. 

p q 2 



All pairs {p,p,q) satisfying (1.151 and (1.161 are called Strichartz pairs. Those for which 
the equality holds in ( |1.16| are called sharp Strichartz pairs. Such estimates first appeared 
in the seminal works |8], |43l HI] and as stated include contributions from, e.g., jTJ, |35] . 
[To] . [26] . and |2l]- 

If one allows variable coefficients, such estimates are well-understood locally-in-time. 
For smooth coefficients, this was first shown in |32] and later for coefficients by |39] 
and gSHMlIlT]. 
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Globally-in-time, the problem is more delicate. Even a small, smooth, compactly sup- 
ported perturbation of the flat metric may refocus a group of rays and produce caustics. 
Thus, constructing a parametrix for incoming rays proves to be quite difficult. At the 
same time, one needs to contend with the possibility of trapped rays at high frequencies 
and with eigenfunctions/resonances at low frequencies. 

Global-in-time estimates were shown for small, long range perturbations of the metric 
in [29] using an outgoing parametrix. In order to keep the parametrix outgoing one 
must allow evolution both forward and backward in time. This construction is based 
on an earlier argument in |48] for the Schrodinger equation. The smallness assumption, 
however, precludes trapping and does not permit a direct application to the current setup. 

On the other hand, a second result of [29] asserts that even for large, long range 
perturbations of the metric one can still establish global-in-time Strichartz estimates 
provided that a strong form of the local energy estimates holds. This switches the burden 
to the question of proving local energy estimates. 

The result in [29] cannot be applied directly to the present problem due to the log- 
arithmic losses in the local energy estimates near the trapped rays. However, it can be 
applied for the near infinity part of the solution. In a bounded spatial region, on the 
other hand, we take advantage of the local energy estimates to localize the problem to 
bounded sets, in which estimates are shown using the local-in-time Strichartz estimates 
of [39], [H]. Thus we obtain 

Theorem 1.4. If (p solves ^g4> = f in Mr then for all nonsharp Strichartz pairs 
(pi,Pi,9i) and (^2,^2, 92) we have 

(1.17) Emj:+)+supE[<P]{d) + ||V0||^.,^-p,„, < i?[0](Sj,) + 11/11' i .,,,4 . 

Here the Sobolev-type spaces H'^'P coincide with the usual H'^'P homogeneous spaces 
in expressed in polar coordinates {r,uj). 

As a corollary of this result one can consider the global solvability question for the 
energy critical semilinear wave equation in the Schwarzschild space. 



(1.18) 



□g0 = ±(f>^ in M 



bo, K(j> = 01 in Eo- 



Theorem 1.5. Let tq > 0. Then there exists e > so that for each initial data ((/jq, (f>i) 
which satisfies 

£;[0](So)<e 



the equation (1.181 admits an unique solution in the region {r > rg} which satisfies the 
bound 

i?[0](E.J-H||0||^.„({,>,„}) <£;[0](So) 
for all indices s,p satisfying 

4 1 1 

- = «+-, < s < -. 
p 2' - 2 

Furthermore, the solution has a Lipschitz dependence on the initial data in the above 
topology. 
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Some further clarification is needed for the function space H^'P{{r > r^}) appear- 
ing above, in view of the ambiguity due to the choice of coordinates. In a compact 
neighbourhood of the center region A4c this is nothing but the classical H^'''' norm. By 
compactness, different choices of coordinates lead to equivalent norms. Consider now the 
upper exterior region Aiji (as well as its three other mirror images). Using the coordi- 
nates {v,x) with X = Lur, we define H'^'P{Aifi) as the restrictions to ]R+ x {|x| > tq} of 
functions in the homogeneous Sobolev space H^'P(R. x K^). 

Acknowledgements: The authors are grateful to M. Dafermos and I. Rodnianski for 
pointing out their novel way of taking advantage of the red shift effect in [H], and to N. 
Burq and M. Zworski for useful conversations concerning the analysis near trapped null 
geodesies. 

2. The Morawetz-type estimate 



In this section, we shall prove Theorem 1.2 We note that the estimate (1.131 is trivial 



over a finite v interval by energy estimates for the wave equation; the difficulty consists 
in proving a global bound in v. By the same token, once we prove (1.131 for some choice 
of To < 2M, we can trivially make the transition to any < 2M due to the local theory. 
Thus in the arguments which follow we reserve the right to take rg sufficiently close to 
2M. 

We consider solutions to the inhomogeneous wave equation on the Schwarzschild man- 
ifold in Mr, which is given by 

Here V represents the metric connection. Associated to this equation is an energy- 
momentum tensor given by 

(9q/3[0] = da(t)df3(j) - ^gapd'^cjydjcj). 

A simple calculation yields the most important property of Qa/s, namely that if solves 
the homogeneous wave equation then Qai3[(f>] is divergence-free: 

More generally, we have 



In order to prove Theorem 1.2 we shall contract Q^p with a vector field X to form 
the momentum density 

Computing the divergence of this vector field, we have 
where 



is the deformation tensor of X. 
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If X is the KiUing vector field K then the above divergence vanishes, 
(2.1) V"Pa[(f),K] = Q if □g(/) = 0. 

This gives rise to the i?o[0] conservation law outside the black hole. 

Naively, one may seek vector fields X so that the quadratic form QapWiT^"^ is positive 
definite. However, this may not always be possible to achieve. Instead we note that it may 
be just as good to have the symbol of this quadratic form positive on the characteristic set 
of Ug. Then it would be possible to make the above quadratic form positive after adding 
a Lagrangian correction term of the form qd'^(f>dj(f>. Such a term can be conveniently 
expressed in divergence form modulo lower order terms. Precisely, for a vector field X, a 
scalar function q and a 1-form m we define 

Pa[(l),X,q,m] = Pa[(f>,X] +q(f)da4'- ^daqcj)'^ + ]^ma<jy^ 

where m allows us to modify the lower order terms in the divergence formula. Then we 
obtain the modified divergence relation 



V"P„[0, X, q, m] = ngcf>(^X(j) + + g[0, X, q, m], 



(2.2) ^ 

, X, q, m] = Q„^[</)]^"^ + qd'^cbd^^P + m^q^d'^cj) + -(V"m„ - V'd^q) 0^ 



Theorem |1.2| is proved by making appropriate choices for X, q and m so that the 
quadratic form Qlcf), X,q,m] defined by the divergence relation is positive definite. In 
what follows we assume that X, q and m are all spherically symmetric and invariant 
with respect to the Killing vector field K. 

Lemma 2.1. There exist smooth, spherically symmetric, K -invariant X, q, and m in 
r > 2M satisfying the following properties: 

(i) X is hounde^, \q{r)\ < r~^, |(7'(?')| < r^^ and m has compact support in r . 

(ii) The quadratic form Q[(j), X,q,m] is positive definite. 



(Hi) X{2M) points toward the black hole, X{dr){2M) < 0, and {m,dr){2M) > 0. 

We postpone the proof of the lemma and use it to conclude the proof of Theorem |1. 2 [ 
Let X, q and m be as in the lemma. We extend them smoothly beyond the event horizon 
preserving the spherical symmetry and the if-invariance. By (2.1 1 we can modify the 
vector field X without changing the quadratic form Q in (2.2 1, 

V°'Pa[(l),X + CK,q,m] = ag(f>(^{X + CK)(j> + q(t)'^ + Q[(l), X, q,m]. 

Here C is a large constant. We integrate this relation in the region 

D ^ {0 < V < vo, r>ro} 



In the (r, v) coordinates 
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using the {r,v,(j) coordinates. This yields 

/ {ng(j) {{X + CK)(I) + qcj)) + Q[(t), X, q, m]) r'^drdvdu 
Jd 



= / {dv, P[(f>, X + CK, q, m])r^drdLU 



v=0 



{dr, P[(t), X + CK, q, m])r^dvdLO. 



We claim that if C is large enough and ro sufficiently close to 2M then the integrals on 
the right have the correct sign, 



(2.3) E[(t)]{vi)<- [ {dv,P[(t>,X + CK,q,m\ydrdio<C £[(!)]{ 

J V=Vi 

(2.4) {dr, P[(^, X + CK, q, m]) > \dr(p\^ + \di(j)\^ + \d^(j)\'' + 0^ 



vi), 



r = ro- 



If these bounds hold then the conclusion of the theorem follows by (ii) and Cauchy- 
Schwarz. 

Indeed, a direct computation yields 
1 



{dv,P[cl),dii]) = 
respectively 



{dr,P[cP,di,]) = \d^cj>\' 



1 - 



2M 



{dr - n' dy)(j)dy(t>. 



On the other hand 

{dv,P[<t>,dr]) =[l-{l- - (^2/x' - (l - d,<l>dr^ 



while 

{dr,P[4>,dr]) 



2m' - (l - |a,</,|2 - (l - ^) \drct^? + 



2Mn 



We compute 

{dv, P[(j), X + CK]) = {X{dv) + C){dv, P[0, dy]) + X{dr){dv, P[<f), dr]). 

For largo onougli C wo have X{dv) + C > C. Thereforo tho first torm on the right is 
negative definite for r > 2M. More precisely, it is only the coefficient of the term 
which degenerates at r = 2M. However, due to condition (iii) in the lemma we have 
X{dr){2M) < 0; therefore we pick up a negative l^^^l^ coefficient at r = 2M. Thus we 
obtain 



-{dv,P[(f),X + CK]) !vC 



+ \dr(pf 



r > 2M. 



Since all the coefficients in the quadratic form on the left are continuous, it follows that 
the above relation extends to r > ro for some ro < 2M depending on C, namely 



(2.5) 



< 2M-ro «: C-^ 
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In order to prove (2.3 1 it remains to estimate the lower order terms P[(j),0,q,m\ in 
terms of the positive contribution above. Since |g| < and m has compact support in 
r, we can bound 

\{dv,P[cP,0,q,m])\ < r-'\cj,\ + ^ + ^-2|^|2_ 

Then by Cauchy-Schwarz it suffices to estimate 



„-2| 1 12„2, 



dr<C- 

which is a routine Hardy-type inequahty. 



/•OO 




2M\ 






[c(, 




+ 1 











\dr(t>rr^dr 



We next turn our attention to \2A\ and begin with the principal part 

(dr, P[(/., X + CK]) = {X{dv) + C){dr, P[(j), d^]) + X{dr)(dr, 9,]). 



Examining the expressions for the two terms above, we see that for rg subject to (2.5 1 
we have 

21VP 
ro 

Next we consider the lower order terms. The contribution of m is 

1 

due to condition (iii) in the Lemma. The contribution of q is 

1 



{dr, P[^, X + CK]) > C\d^<j>\^ + \d^cb\^ - (l - ^) r = tq. 

contribi 
-(m, dr)(f>'^ > (jP' 



q(j){dr,d(j)) ~ ^(f^ {dr , dq) . 
The coefficient of the second term is 

(l - which is negligible for rg close to 2M. 

In the first term we have 

(dr,d(/)) = (i - (i - (i - 

All terms involving ^1 — are negligible, and since q is bounded we get 
for large enough C. 



Proof of Lemma \2.1\ It is convenient to look for X in the {r,t) coordinates, where we 
choose the vector field X of the form 

X = Xi+5X2, 5<1 

where 

and a and 6^1 — will be chosen to be smooth. Note that X is a smooth vector field 
in the nonsingular coordinates (r, w), since in these coordinates we have 



X, = a{r) n 1 - + q\ , = b{r){l ^) 



dr 
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We remark that the vector field X2 is closely related to the vector field Y introduced 
earlier in in order to take advantage of the red shift effect. However, in their con- 
struction y is in a form which is nonsmooth near the event horizon and which is restricted 
to the exterior region. 

The primary role played by X2 here is to ensure that X + CK is time-like near the 
event horizon. The red-shift effect largely takes care of the rest. 

For convenience, we set 

2M\ 1 



'.M = (i-i^i)i*(Aw) 



A direct computation yields 
(2.6) 

respectively 



2M^ 



a'{T){dr(jyf + a{r)- 



3M 



-Mr)d^4>d^4> + X^<t>Og^, 



\yir) ((1 



2M\ 



]dr(j) ~ dt(f) 



(2.7) 



where 



2A/X , 



r / 



6'(r) 1)701 



1 



2M 
r J 



(dt<i>f 



1 



2M 



We choose a so that the first line of the right side of (2.61 is positive. This requires that 
(2.8) 



a'(r) 



> 



a(3M) 0. 



We choose h so that the first line of the right hand side of (2.7 1 is positive. Precisely, we 
take h supported in r < 3Af with 



2M 



1 



r G [2Af,3M]. 



with bo smooth, decreasing in [2M, 3M) and supported in {r < 3M}. In particular 
that guarantees that bo{2M) > 0, which is later used to verify the condition (iii) in the 
Lemma. 

The exact choice of bo is not important, and in effect b only plays a role very close 
to the event horizon r = 2M. Even though b is singular at 2M, the second term of the 
coefficient of I'^'Pl'^ in the second line of (2.7 1 is nonsingular. Hence if 6 is sufficiently 
small this term is controlled by the first line in (2.6 1. 
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Taking the above choices into account, we have 

r 



(2.9) 



where 



1 



O 



) a'{r){dAY + 5\b'{r){^[ 
(r - 3M)2 



2M\ 
r J 



19^0 - dt(t) 



2M\ 
r / 



The last term in (2.91 is a Lagrangian expression and is accounted for via the q term. 
The first three terms give a nonnegative quadratic form in V0. This form is in effect 
positive definite for r < 3M, where b' > 0. However for larger r it controls dr4> and ^0 
but not 9i0. This can be easily remedied with the Lagrangian term. Precisely, we choose 
q of the form 

(r-3M)2 

q = qo + diqi, qi[r) ^ X{r>5M/2} • 

where X{r>5M/2} is a smooth nonnegative cutoff which is supported in {r > 5M/2} and 
equals 1 for r > 3M. The positive parameter Si is chosen so that Si ^ S. Then the only 
nonnegligible contribution of Siqi is the one involving dt4>. We obtain 



(2.10) 



2M 



J' 



S-b'ir) 



2M\ 
r J 



+0 



(r- 3Af)^ 



- 9f0 



The contribution of qi can be made arbitrarily small by taking S\ small. Hence it will 
be neglected in the sequel. At this stage it would be convenient to be able to choose o 
so that V"(9aii(r) < 0. A direct computation yields 



V"aaii(r) = -La 



with 



M.) = ->..[(! -^).'*{(i-^)5i.*(.M.-))} 



Unfortunately it turns out that the condition La > and (2.81 are incompatible, in the 
sense that there is no smooth a which satisfies both. However, one can find a with a 
logarithmic blow-up at 2M which satisfies both requirements. Such an example is 

-2My 



a(r) = r-^ (r - 3A//)(r + 2M) + GA/^ log 



M 



This is in no way unique, it is merely the simplest we were able to produce. One verifies 
directly that 

a'(r-) > La(r) > r^^ 

To eliminate the singularity of a above we replace it by 

a,(r) = \mR) 
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where e is a small parameter, 

/ r — 2M 

R^{r~ 3A/)(r + 2M) + log 

and 



M 



where / is a smooth nondecreasing function such that f{R) — R on [—1, oo] and / — —2 
on (— oo, —3]. The condition (2.81 is satisfied uniformly with respect to small e; therefore 
the choice of 6 is independent of the choice of e. 

With this modification of a we recompute 

La, = neR)La + 0{e)f"{eR) + O ( f 1 - — ) ') f"{eR). 



This is still positive except for the region {eR < —1}. To control it we introduce an m 
term in the divergence relation as follows. 

Let 7(r) be a function to be chosen later. We set 

2 



Then 



while 



mt = 6b'{r)(^l — ^ 7, mr = Sb'{r)(^l ^ 0. 

m^d^cj, = 56'(r)7(r) (l - ^) ({l - ^)dr(b - dtA , 



Hence, completing the square we obtain 

Q[<l>,X,q,m] = (1 - g^)V'^)(9.0r + o( ^''^f^^' ) m 

+ Sl^ir) (^(1 - - d,cl>+ (1 - ^)7^ 

where the coefficient n is given by 

1, 0/ 2M\„ / , , 2M\\ M(r) / 2M\^ , 

We assume that 7 is supported in {r < 3M} and satisfies 

0<7<1, 7' > -1. 

Then for r > 3M we have 
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while for r < 3M we can write 

La, + Sy{r)(l~^^\'{r) + 0{S). 

K ei? > — 1 then, using the bound from below on 7', we further have 

n>La + 0{S) 

which is positive provided that S is sufficiently small. On the other hand in the region 
{eR < —1}, we have 

n > ls[l - '^J^)\'{r)i{r) + 0{5) + 0{e)r{eR) + O (^e^[l - ^) /'"(ei?). 

The 7' term can be taken positive, while all the other terms may be negative so they 
must be controlled by it. The restriction we face in the choice of 7' comes from the fact 
that < 7 < 1. Hence we need to verify that 

/ S + e\neR)\ + e'(l-—y\r'{eR)\«S. 

JeR<-l ^ r J 

Indeed, the interval of integration has size < e"'^^ ^; therefore the above integral can be 
bounded by 

which suffices provided that e is small enough. 
Finally, note that 

X{dr){2M) = (^a{r)(l - ^) +,56(r)(l - (2Af) < 0, 

{m,dr){2M) = (^Sb'{r)(^l - — ) 7^ {2M) > 0. 
So (iii) is also satisfied. □ 

3. Log-loss local energy estimates 

The aim of this section is to prove a local energy estimate for solutions to the wave 
equation on the Schwarzschild space which is stronger than the one in Theorem |1.2[ 
Consequently, we strengthen the norm LEq to a norm LE and we relax the norm LEq 
to a norm LE* which satisfy the following natural bounds: 

(3.11) ll'/'llU,<ll'/'lliB<ll'^llU„ 

respectively 

(3.12) \\f\\lEl,<\\f\\lE'<\\f\\lES- 

We note that these bounds uniquely determine the topology of the LE and LE* spaces 
away from the photon sphere and from infinity. This is due to the fact that the local 
energy estimates in Theorem |1.2| have no loss in any bounded region away from the photon 
sphere. To define the LE, respectively LE* , norms we consider a smooth partition of 
unity 

1 = 



= Xeh{r) + Xps{r) +Xoo{r) 



STRICHARTZ ESTIMATES ON SCHWARZSCHILD BLACK HOLE BACKGROUNDS 



19 



where Xeh is supported in {r < 11A//4}, Xps is supported in {5M/2 < r < 5M} and Xo 
is supported in {r > 4M}. Then we set 

(3.13) UWIe = WXehmE,, + WXps^Ie,^ + llXoo^lliiJ,,, 

respectively 

(3-14) = IIXeh^llLj, + IIXp^^llL. + \\Xoo(t>\\lE:,- 



The norms LEps and LE*^ near the photon sphere are defined in Section 



3.1 



(3.201, respectively (3.211; their topologies coincide with LEm, respective^ 



below, see 
YtElj, away 



from the photon sphere. 



With these notations, the main result of this section can be phrased in a manner 
similar to Theorem 11.21 



Theorem 3.2. For all functions cj) which solve Ogcj) = f in Mr we have 
(3.15) snpEmd) + Em^+) + < i?M(S^) + UWIe-^. 

v>0 

We continue with the setup and estimates near the photon sphere in Section |3.1[ the 
setup and estimates near infinity in Section |3.2| and finally the proof of the theorem in 
Section |331 



3.1. The analysis near the photon sphere. Here it is convenient to work in the 
Regge- Wheeler coordinates given by 

r* r + 2M log(r - 2M) - 3M - 2M log M. 

Then r — 3M corresponds to r* — 0, and a neighbourhood of r = 3M away from infinity 
and the event horizon corresponds to a compact set in r* . In these coordinates the 
operator Dg has the form 

(3.16) r(^l-^^ n.r-i - Lrw = df - d% - "^^^d^ + V{r), V{r) = r-'d%r. 

For r* in a compact set the energy has the form 

« y (9*0)2 + (5,. 0)2 + {d^cj>ydrdLu, 
and the initial local smoothing norms are expressed as 



respectively 



On the other hand 



l/llifiS « / r*-'fdrdu;dt. 



12 ^ I fa J,^2 I /a J,^2 I fa J,^2 , j,2 

\lEm ' 



leIj ~ J fdrdujdt. 



20 JEREMY MARZUOLA, JASON METCALFE, DANIEL TATARU, AND MIHAI TOHANEANU 



In the sequel we work with spatial spherically symmetric pseudodifferential operators 
in the {r*,uj) coordinates where u; G S^. We denote by ^ the dual variable to r* , and by 
A the spectral parameter for {—di^) s . Thus the role of the Fourier variable is played by 
the pair A), and all our symbols are of the form 

To such a symbol we associate the corresponding Weyl operator A^. Since there is no 
symbol dependence on iu\ one can view this operator as a combination of a one dimensional 
Weyl operator and the spectral projectors 11;^ associated to the operator (— As2)2j namely 

X 

All of our estimates admit orthogonal decompositions with respect to spherical har- 
monics, therefore in order to prove them it suffices to work with the fixed A operators 
a™ (A), and treat A as a parameter. However, in the proof of the Stricliartz estimates 
later on we need kernel bounds for operators of the form A™, which is why we think of 
A as a second Fourier variable and track the symbol regularity with respect to A as well. 
Of course, this is meaningless for A in a compact set; only the asymptotic behavior as 
A ^ 00 is relevant. 



1 y<i, 
y y>2. 



Let 7o : M — > M"*" be a smooth increasing function so that 

lo{y) 

Let 7i : M+ — > M+ be a smooth increasing function so that 

71(2/) = 



y2 y< 1/2, 
1 y> 1. 

Let 7 : — > 1R+ be a smooth function with the following properties: 

r 1 z<c, 

7{y,z)=l 70(2/) y<^/z/2,z>C, 
[ zhiiy^z) y > z > C, 

where C is a large constant. In the sequel 2; is a discrete parameter, so the lack of 
smoothness at 2; = C is of no consequence. 

Consider the symbol 

aps{r*,^,X) = 7(-ln(r*2 + A-V),liiA), 

and its inverse 



' 7(-ln(r*2 + A-2^2)^inA)- 

We note that if A is small then they both equal 1, while if A is large then they satisfy the 
bounds 

^^^^^ l<aps(r*,?,A) <ap,(r*,0,A) < (InA)^, 

(lnA)-5 < a-\r\0,X) < a-\r\^,X) < 1. 
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We also observe that the region where > z corresponds to r*^ + A^^^^ < e^^'"'^. 
Thus differentiating the two symbols we obtain the following bounds 

(3.18) |9,".9f9>p,(r*,^,A)| < c„,/3,,A-'3"''(r* V A'^f + f-^)-^ , 
respectively 

(3.19) |a,".9fa>-i(r*,C,A)| < c^,p,,a-^{T\^,\)\-P-'^{r*^ + \-^e + e-^^r^ , 

where a + [3 + v > Q. These show that we have a good operator calculus for the cor- 
responding pseudodifferential operators. In particular in terms of the classical symbol 
classes we have 

o-ps,ap} ^ SIq, (5 > 0. 



Then we introduce the Weyl operators 

A„., = ^a»(A)nA, 

A 



respectively 



By (3.181 and (3.191 one easily sees that these operators are approximate inverses. More 
precisely for small A, In A < C, they are both the identity, while for large A 

ll«p.(A)(a;/)'^(A) - /lU.^i. < \-^e^\ InA > C. 
Choosing C large enough we insure that the bound above is always much smaller than 1. 
We use these two operators in order to define the improved local smoothing norms 

(3.20) uue,^ = \\A-}nHi 

(3.21) WfWLE' = I 



\A-}Vt,A\\L^ 



\Apsf\\ 



L2. 



Due to the inequalities (3.171 we have a bound from above for a™ (A) 



|a;s(A)/|| 



L2 



< 



||ap,(r*,0,A)/||i2 < |||ln|r* 



ps 

L2, 



respectively a bound from below for (flp^ )"'(A) 



ll«r(A)/|U^ > ||a;,i(r*,0,A)/|U. > |||ln|r*rVl| 



L2 



for / supported near r* = 0. In particular this shows that for / supported near the 
photon sphere we have 



(3.22) 



m 



> II link* 



-1V0II 



L2, 



l/ll 



LE' 



< 



|ln|r1|/|U. 



which makes Theorem |1.3| a direct consequence of Theorem |3.2[ 
Our main estimate near the photon sphere is 

Proposition 3.3. a) Let <j) he a function supported in {5M/2 < r < 5M} which solves 
□ = /. Then 



(3.23) 



\le^ 



<ll/l" 



LE' 
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b) Let f e LE*^ be supported in {llM/i < 7' < 4Af}. Then there is a function (p 
supported in {5M/2 < r < 5M} so that 

(3.24) supSM + ll^llL^^ + llDg^-ZllL. <|1/|1L.^. 



Proof. Due to (3.161 we can recast the problem in the Regge- Wheeler coordinates. De- 
noting u — r(j), g = ^1 — ^^r/ we have Lj^wu — 9- Also it is easy to verify that for 
and / supported in a fixed compact set in r* we have 

UWle^, ~ \\u\\le^,, II/IUs;, ~ hhE;^- 
Hence in the proposition we can replace and / by m and g, and Dg by Ljiw 

To prove part (a) we expand in spherical harmonics with respect to the angular variable 
and take a time Fourier transform. We are left with the ordinary differential equation 

(3.25) {d^^,+Vx^r{r*))u = g, 

where 



5 \ + V. 



Depending on the relative sizes of A and r we consider several cases. In the easier 



cases it suffices to replace the bound (3.231 with a simpler bound 

(3.26) \\drM\L- + i\T\ + \\\)\\u\\L. 



< 



Case 1: A,r < 1. Then we solve (3.251 as a Cauchy problem with data on one side 
and obtain a pointwise bound , 

\u\+\Ur* \ < 



which easily implies (3.261. 



Case 2: A ^ r. Then V\,rir*) w r for r* in a compact set; therefore (3.251 is 
hyperbolic in nature. Hence we can solve (3.251 as a Cauchy problem with data on one 
side and obtain 

t\u\ + \Ur' \ < \\g\\L^, 



which implies (3.261. 



Case 3: A r. Then Vx^rir*) ~ ~A^ for r* in a compact set; therefore (3.251 is 
elliptic. Then we solve (3.251 as an elliptic problem with Dirichlet boundary conditions 
on a compact interval and obtain 

A5|u|+A5|u,.| < \\g\\L2, 



which again gives (3.261. 

Case 4: A ~ r ^ 1. In this case (3.261 is no longer true, and we need to prove (3.231 
which in this case can be written in the form 

(3.27) Wdr^uh^ + XUa-lr{X)uh. < ||a- (A)5||l=, 

where u,g are subject to (3.251. The dr'U term above is present in order to estimate the 
high frequencies |^| ^ A. For lower frequencies it is controlled by the second term on the 
left of K27\. 
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The potential V in ( |3.25[ l can be treated perturbatively in (3.231 and is negligible. The 
remaining part of V\^r{r*) has a nondegenerate minimum at r = 3M which corresponds 
to r* = 0. Hence we express it in the form 

where W is smooth and has a nondegenerate zero minimum at r* = and |e| < 1. 
We now prove the following: 

Proposition 3.4. Let W be a smooth function satisfying W{0) = W'{0) = 0, W"{0) > 0, 
and \e\ < 1. Let w be a solution of the ordinary differential equation 

d% + \^{W{r*) + e))w{r*) ^ g 



supported near r* = 0. Then (3.27 1 holds. 



It would be convenient to replace the norm on the right in (3.271 by || aps(r*, 0, X)g\\]^2. 
This is not entirely possible since this is a stronger norm. However, we can split g into a 
component gi with aps{r* ,0, X)gi G L^ plus a high frequency part: 

Lemma 3.5. Each function g E supported near the photon sphere can be expressed 
in the form 

9 = 91 + X^^d^,g2 
with gi and 52 supported near the photon sphere so that 



(3.28) ||ap,(r*,0,A)5i|U2 + |||r*- + e 



-\/ln Ai 



\92||l2+A ^119^-3211 



L2 



< 



I<s(A)5||l- 



Proof. The symbols aps{r* ,0, X) and aps(r*,^, A) are comparable provided that 

ln(r*^ + e-^^) « ln(r*^ + e"^^ + A^^^^)^ 
This includes a region of the form 

D = (ln(A-2f ) < 1 ln(r*' + e"^) 



We note that the factor |, arising also in the exponent of the second term in (3.281, is 
somewhat arbitrary. A small choice leads to a better bound in (3.281. 

If X is a smooth function which is 1 in (— cx),— 1] and in [0,oo) then we define a 
smooth characteristic function x_d of the domain D by 

1 



XD(r*,C,A) = x(ln(A"2e')--ln(r* 
One can directly compute the regularity of xd, 

XD e Sis, S > 0. 

To obtain the decomposition of g we set 

92 = q'^g, 

where the symbol of q is 

q{r*,^,X) = X'r^{l-XD). 



/In A 



))■ 
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Since (cips)^ is an approximate inverse for a™^, the estimate for 52 in the lemma can be 
written in the form 

(3.29) \\{r*' + e-^^)5g»(a-i)-(A)/|U. + A-^|19.-g-(a-i)-(A)/|U. < ||/|U.. 

In the first term it suffices to look at the principal symbol of the operator product since 
the remainder belongs to OPS^l'^^ for all S > 0. To verify that the product of the 
symbols is bounded we note that a~/ is bounded. For the other two factors we consider 
two cases. If |^| > A then both factors are bounded. On the other hand if |^| < A then 
in the support of q we have 

which gives 

9< (^*' + e-^)-5. 



The estimate for the second term in (3.291 is similar but simpler 



It remains to consider the bound for gi, which is given by 

I 

As above, the bound for gi can be written in the form 

\\a,,{r*,Q,X){l + X-^Dl,q-){a^}r{\)f\\L. < 

The three operators above belong respectively to S( g, g, and S'f ^ for all 6 > 0. Hence 
the product belongs to S( g, and it suffices to show that its principal symbol is bounded. 
But the principal symbol of the product is given by 

aps{r*,0,X)xDa~s{r*,tX), 

which is bounded due to the choice of D. 

Finally we remark that as constructed the functions gi and 172 are not necessarily sup- 
ported near the photon sphere. This is easily rectified by replacing them with truncated 
versions, 

51 := Xi(''*)5i, 92-^x{r*)g2 
where xi is a smooth compactly supported cutoff which equals 1 in the support of g. It 



is clear that the bound (3.281 is still valid after truncation. 



□ 



Using the above decomposition of g we write u in the form 

U = A^^(72 + u. 

For the first term we use the above lemma to estimate 

A||A-2^2||l= + \\X''dr*g2\\L- < ||a;,(A)g|U2, 
which is stronger than what we need. For u we write the equation 
(3.30) + X'^iW + e))u = 5, g = .91 - (VK + e).g2. 

For g we only use a weighted bound, 

\\{X-' + \W + e\)-i~gU.<\\a-AX)g\\L^, 
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which is obtained from the weighted bounds on gi and 52 by Cauchy-Schwarz. 

For u on the other hand, it suffices to obtain a pointwise bound: 
Lemma 3.6. For each < a < 1 and each function u with compact support, we have 
A||(%/r(A)w|U2 < \\{a+\W + e\)-'-dr,u\\L«'+X\\{a+\W + e\)--u\\Loo. 

Proof. Since W has a nondegenerate zero minimum at 0, if e > —a then a + \W + e| « 
cr+lel + W. Hence without any restriction in generality we can replace (e, cr) by (0,(7+ |e|). 
Thus in the sequel we can assume that either e = or e < —a. We consider three cases: 

Case 1: \e\,a < e^^'"'^. We consider an almost orthogonal partition of u in dyadic 
regions with respect to r*: 

So<S<l 

For each piece we can freeze r* in the symbol of o~/ and estimate in 

\\{a;}r{X)u,h'^ « \\a;}{s,D,\)U.s\\L^, || (a^i)"'(A)zi<.JU2 « ||api(0, Z), A)u<,JU2. 

In the first case we use the symbol bound 

Aa-'(s,C,A)< A|lns|-i + s-^|e|, 5 > 0, 

where the second term accounts for the region where |^| > Xs^ . This yields 

A||(a-i)«'(A)w3|U2 < \lns\-^X\\us\\L-+s-'\\dr,Us\\L-. 

In the support of Ug we have cr + \W + e| w s^; therefore by Cauchy-Schwarz we obtain 

A||(a,-;)-(A)«,|U. < I lns|-^A||(a + \W + e\)iu\\L^ + s'-'\\{a + \W + e|)-^5..u|U^ . 

The summation with respect to s follows due to the almost orthogonality of the 
functions {a~^)^{X)us, 



X 




< X\\{a +\W + e|)iu||ioo + \\{a+\W+ e\)-Wr^u\\Lo 
This orthogonality is due to the fact that the kernel of {a~g)^{X) decays rapidly on the 



X scale in r* , therefore the overlapping of the two functions of the form [a^g)"^ {X)us 



is trivially small for nonconsecutive values of s. 

On the other hand for the center piece tt<so a similar computation yields 



A||(a;;)-(A)S<,JU. < \ lnso\-'X\\u^so\\L^+So'\\dr'U<so\\ 



L2 



< X\\{a+\W + e\)-^u\\L^+sl-'\\{a+\W + e\y-^dr^u\\L^ 



where the weaker bound in the first term is due to the fact that 

{a+\W + e\)-^dr* <lnA. 



|r*|<so 
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<^ase 2: e = 0, cr > e"^'"^. Then we select sq by Sq = ct and partition u into 

Sq<S<- 1 

The analysis proceeds as in the first case, with the simplification that there is no longer 
a singularity in the weig ht {a+\W + e\)-2 for \r*\ < sq. 

Case 3: a, e"^^^^ < — e. Then we select so by Sq = — e and partition u into 



'0 



U = U<so + "so + 

sa<s<l 

Then all pieces are estimated as in Case 2 with the exception of Usg, where we have to 
contend with the singularity in the weight. However, compared to Case 1 we have a 
better integral bound 

/ \W + e\-^dr* <l 
and the conclusion follows again. □ 

Due to the above lemma, it suffices to prove that, given |e| < 1 and {u,g) supported 
near the photon sphere so that 

(3.31) {d^,+X^{W + e))u^g, 
then for some A^^ < cr < 1 we have 

(3.32) \\{a + \W + e\)-idr'u\\L^ + X\\{a + \W + e\)iu\\L^ < \\{X-^ + \W + e\)- g\\ li . 



We remark that the first term on the left gives the bound for Ur* in (3.271. 
We consider three subcases depending on the choice of e: 

Case 4 (i) e A^^. Then \W + e\ ~ r*^ + e. Choosing cr = e, it suffices to prove 
that: 



Lemma 3.7. Suppose that e 3> A ^. Then for u with compact support solving (3.311, 
we have 

A(e + r*^)3|w| + {e + r*'^)-^\ur* \ < ||(e + r*^)"i.g||ii. 

Proof. We solve ( 3.31[ l as a Cauchy problem from both sides toward 0. For this we use 
an energy functional which is inspired by the classical WKB approximation. 



E{u{r*)) = A^(Ty + e)5w^ + {W + e)~^-ui, + -Wr-{W + ey 



UUr 



2 

Since \Wr-' \ ^ W^, the condition e ^ A~^ guarantees that E is positive definite. Com- 
puting its derivative, we have 

^Eiuir*)) =i(Ty,.,.(W^ + e)-i - ^W^,(W + e)-l)uur. 



This leads to the bound 



+ {2{W + ey^Ur' + ^Wr^iW + e)'^u)g. 



< X-^{W + e)-^ E{u{r*)) + E^u{r*)){W + e)-i g. 
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Since 

the conclusion follows by Gronwall's inequality. 

Case 4 (ii) |e| < A^^ We choose cr = A"^ Then a H 
suffices to prove the pointwise bound: 



\W- 



e\ w A-i +r* 



□ 



Hence it 



Lemma 3.8. Suppose that |e| ^ A ^. Then for u with compact support solving (3.3f I, 
we have 

A(A-i + r*Y^\u\ + (A-i + r*Y^Ur'\< \\{X-' + r*')-'^ g^^. 

Proof. We use the same energy functional E as above. However, this time E is only 
positive definite when W ^ A~^ or equivalently \r*\ ^ A~2. Applying Gronwall as in 
the first case yields the conclusion of the Lemma for |r*| ^ A~2 . 



In the remaining interval {I?"*! ^ A 2 } we view (3.31 1 as a small perturbation of the 
equation d^,u = g. Precisely, we can use the energy functional 



£;iKr*)) = Ai|u|2 + A5|w,.|2, 



which satisfies 



^-E,{u{r*)) < \--E,{u{r*)) + Ef {u{r*))\-^ g. 
This allows us to use Gronwall's inequality to estimate the remaining part of u. 



□ 



Case 4 (iii) — e 3> A ^. Then we choose cr = |e| 3 A 3 and prove the pointwise bound: 



Lemma 3.9. Suppose that — e 3> A ^ . Then for u with compact support solving (3.31 1 
we have 

\{\W + e\ + \e\^\-i)-^\u\ + {\W + e\ + \e\^\--^)--^\urA< \\{\W + e\ + \€\^\~i)--^g\\Li. 

Proof. The energy E above is still useful for as long as it stays positive definite, i.e. if 
(3.33) \Wr.\<2X\W + e\^l'^. 

Given the quadratic behavior of W at 0, this amounts to 

+ e > A"3|e|3. 



In this range, due to (3.331 we obtain, as in Case 4(i), 
d 



dr 



■E{u{r*)) < X-'^iW + e)-^E{u{r*)) + E^{u{r*)){W + e)-i< 



which by Gronwall's inequality and Cauchy-Schwarz yields the desired bound. 
In a symmetric region around the zeroes oi W + e, 

\W + £\ < A-s|e|s, 



the bounds for u and Ur- remain unchanged and the equation (3.31 1 behaves like a small 
perturbation of d^^u = g, and we can use a straightforward modification of the argument 
above. 
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Finally, in the region 

[r_,r+] ={W + e< -CX-ie^}, 
we use an elliptic estimate. Denote 



uj=\W + e\ + |e|5A~3. 



Then multiplying the equation (|3.31|l by —Xu and integrating by parts we obtain 
(3.34) 



X\dr-u\'^ + X^\W + e\\u\'^dr* = 



-Xugdr* + Xuu^ 



< 



X\\ujiu\\L'^(^r^.r+)\\^^ ^ g\\ + \\^^ * gWh , 



where for the boundary terms at r± we have used the previously obtained pointwise 
bounds. On the other hand from the fundamental theorem of calculus one obtains 



< 



Xldr-ul"^ + X^\W + e\\u\^dr* 



where the bound (3.331 is used for the derivative of W in [r_,r_|_]. Combining the last 
two inequalities gives the desired bound for u, 



(3.35) A| 
Returning to (3.34 1, it also follows that 



(3.36) 



Xldr'ul'' + X'^IW + e\\u\^dr* < 



It remains to obtain the pointwise bound for Ur*. In [r„, r+] we have W < \e\ therefore 
^ Given Tq G [r_,r_|_] we consider an interval G / C [r_,r+] of size 

|/| — cX~^uj{rQ)~2 with a small c. In / the size of the weight oj is constant; indeed, oj 
can change at most by 

\I\\Wr*\ = cX-^uj{r*)~h^ < cuj{r*) 
where at the last step we have used the bound a;(r5) > jej^A^i. 



Within / we first use the bound (3.361 to estimate the average u^, of u^* in /, 

\uU' <\ir' JWr-fdr* <u;ir*)^\\c.-ig\\l 

It remains to compute the variation of u^* in I, which is estimated using the equation 
d^, + X^{W + e) = 5 and ( [slsf , 

\d^,,u\dr* < X^uj\u\ + \g\dr* < a>(r*)3 H^-i^H^i 

Together, the last two bounds show that 

\ur'{r^)\ < ^{r*o)i\\i^-i9\\L^- 
The proof of the lemma is concluded. □ 
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We continue with part (b) of the proposition. We switch to the Regge- Wheeler co- 
ordinates. By taking a spherical harmonics expansion it suffices to prove the result 
at a fixed spherical frequency A. Let g\ be at spherical frequency A with support in 
{llAf/4 < r < 4M}. Using a time frequency multiplier with smooth symbol we can 
split g\ into two components, one with high A) time frequency and one with low time 
frequency. We consider the two cases separately. 

Case I. g\ is localized at time frequencies {|r| ^ (1 + A)}. This corresponds to Cases 
1,2,3 in the proof of part (a). As a consequence of the results there we have the a-priori 
bound 

i\r\ + X)\\u\\L^<Ud-^,+V^,M\L^ 
for all u with support in {5M/2 < r < 5M}. By duality this implies that for each g & 
with support in {5A//2 < r < 5M} there exists a solution ti to 

(0% + Vx,r)v = g 

iif]{5M/2 < r < 5M} with 

(M + A)||HU. <|l5||L- 
Applying this at all time frequencies |r| 3> (1 + A) we find a solution u\ to 

(3.37) LrwUx = gx 

in {5M /2 < r < 5M} so that 

{l + \)\\ux\\L- + \\dtUx\\L-<\\gx\\L^- 



Multiplying equation (3.371 by Xps'^>^ ^^'^ integrating by parts we obtain 

\\dr'{XpsU\)\\l2 < A^llxpsMAllis + \\XpsdtUx\\l2 + \\ux\\l2 + hxWh- 
Hence the function vx — XpsU\ satisfies 
(3.38) W^vxWl- < \\9x\\l-- 

On the other hand, since gx is supported in the smaller interval {llAf/4 < r < 4M}, it 
follows that Vx solves the equation 

LrwV\ - g\ = [Lrw,Xps]ux- 

Here the right hand side is supported in a region, away from the photon sphere, where 
the and LE*^ norms are equivalent. Then this is seen to satisfy 

ll^w^A - gxh^ < ||5a||l2 



by applying (3.381 with Xps replaced by a cutoff with slightly larger support. 

Finally, the standard energy estimates for vx allow us to obtain uniform energy bounds 
for Vx from the averaged energy bounds in (3.381, thus improving (3.381 to 

(3.39) IIVi-aIIl^ + ||Vwa||l~l2 < ||.9a||l2. 

Case II. gx is localized at time frequencies {|t| ^ (1 + A)}. This corresponds to Case 
4 in the proof of part (a). We first observe that the result in part (a) can be strengthened 
to 

(3.40) II0IIW. < ll/llL;^+LiL- 
■'no boundary condition is imposed on v 
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Indeed, suppose that f — fi + f2 with fi e LE*^ and /2 G L^L^. We solve the forward 
problem 

By Theorem |1 .2 1 and Duhamel's formula we have 



||02||l£o ^ II/2||lil2- 

We truncate (f>2 Xps(?')</'2 in a slightly larger set than the support of /2 and compute 

W^giXpsh)- f2\\LE;^, = ||[n3,Xps]'/'2||L£;;^ « IPg, Xps]02 1| < |1 <?!>2 || LBo 

since the above commutator is supported in a compact set in r away from the photon 
sphere. 

From Duhamel's formula and part (a) of the proposition it follows that 

\\Xps<p2\\LEp, < II/2||l1L2- 

On the other hand applying directly part (a) of the proposition to — Xps'p2 we obtain 

110- Xps02||L£;p, < \\Og{(l)- Xps(I>2)\\le;^ < II/iIIlb;, + II/2IU1L2. 
Hence ( |3.40| follows. 

As a consequence of ( |3.40 1 we obtain 

A||(a;/r(AK|U. < inf (||a;,(A)5i|U2 + II52IILIL2) . 

By duality, from this bound from below for Lfiw, we obtain a local solvability result. 
Precisely, for each g\ at spherical frequency A with support in {5M/2 < r < 5M} there 
is a function u\ in the same set which solves 

(3.41) Lrwux = 5a 
and satisfies the bound 

(3.42) A(||(a-i)-(A)u,|U. + ^aIIl-lO < \\a;M)gx\\L^. 

Since (a~/)™ has an inverse in OPS^q, from the first term above we also obtain an 
bound for u\, namely 

(3.43) ^'-'\\ux\\L^<\\a;,i\)gx\\L^. 

Since gx is localized at time frequencies |t| < (1 + A), it follows that ux above can be 
assumed to have a similar time frequency localization. Hence (3.421 also gives 

(3.44) \\ia;^'riX)uxt\\L^ + Wuxth-L^ < \\a^,{X)gx\\L^ . 

We can also obtain a similar bound for the r* derivative of ux- For the local energy part 
we multiply (3.411 by Xps{{0'ps)"' {^))'^XpsUx- After some commutations where all errors 
are bounded using the previous estimates we obtain 

\\{ap}r{X)dr.{xpsUx)\\h< X'\\{a;}r{\)xpsux\\h + Wi^-^r WXpsUxtWh 

+ >^'-''\\ux\\l. + \\9x\\h. 

For the L°°L'^ bound on dr-'{xpsUx) we consider a smooth compactly supported function 
x{t)- Then multiplying ( |3.41| by X^Xp^UA and commuting we obtain 



\\xdr'{XpsUx)\\\2 < A^||xXps"A||i2 + ||xXpsWAt|li2 + ||uA||i2 + llffAll 



L2- 
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Taking also (3.421 and (3.44 1 into account we have a bound on local averaged energy for 
XXpsUy. 



\\dr'{xXpsUx)\\l2 + X'^\\xXpsU\\\l2 + WdtixXpsUx] 



K(A)5A||i.. 

By energy estimates applied to xXps^x we can convert the averaged energy bound into 
a pointwise energy bound to obtain 

\\dr'{xXpsUx)\\l^L2 + X^WxXpsUxWl^L^ + \\dt{xXpsUx)\\l^L^ < \\ap,{\)gx\\l2. 

Summing up ( 3.42[ |, (3.44 1 and the similar bounds above for the r* derivatives we finally 
obtain 

\\{a;}riX)\/{xpsUx)\\L^ + \\V{xpsUx)\\l-l^ < \\9x\\le', 
where V — {dr* ,dt,X). 

On the other hand if gx is supported in {llAf/4 < ?- < 4Af} then ux solves the 
equation 

LrwXpsUx — 9x = [Lbw,Xps]ux- 
The right hand side is supported away from the photon sphere, where the and LE^ 
norms are equivalent. Then, by applying Theorem 
slightly larger support, this is seen to satisfy 

\\LrwXpsUx- gxllLE^ < \\g\\\LE;^- 
The proof of the proposition is concluded. □ 



1.2 



with Xps replaced by a cutoff with 



3.2. The analysis at infinity. In the Schwarzschild space A4, if a function u in is 
supported in {r > AM} we interpret it as a function in KxM"^ by setting u{t, x) — u{t, r, ui) 
for X — ru. We now state the analogue of Proposition |3.3| 



Proposition 3.10. a) Let (p solve Ogcj) = in {r > 4M}. The 



IIXoo'/>| 



2 

LEn 



< 



2 

\lEo 



b) Let f G LE^j be supported in {r > 4M}. Then there is a function (j) supported ir. 
{r > 3M} which solves Dgcj) — f in {r 3> M} so that 



(3.45) 



llE.s + W^s'l^-fWh^WfWlEt 



Proof, a) For i? > we denote by x>R ^ smooth cutoff function which is supported in 
{|a;| > R} and equals 1 in {\x\ > 2R}. If i? > 4M then 

IKXco - X>R)(t>\\lEM - ll'?^lliiSo- 

It remains to show that for a fixed sufficiently large R we have 

\\x>RmE.,<\miEo + E[4'm. 

For this we notice that x>r4' solves the equation 

(3.46) Og{x>R(t>) = /i(2:)V0 + /2(x)0, 

where /i and /2 are supported in {R < \x\ < 2R}. If R is sufficiently large then 
outside the ball {|a;| < R} the operator Dg is a small long range perturbation of the 
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d'Alembertian. Then the estimate (1.5 I appUes, see e.g. pOi Proposition 2.2] or |28l 
(2.23)] (with no obstacle, = 0) and we have 



\X>BmLE 



' <E[x>Rcf>m + \\09ix>BM 



<Emo) + \\[ag,x>Rl 

< mm 4 



2 

LBn' 



where in the last two steps we have used the compact support of Og(x>R4>) = i^g^ X>r]4>- 

b) Let R be large enough, as in part (a). For > R the Schwarzchild metric g is 
a small long range perturbation of the Minkowski metric, according to the definition in 
|29) . We consider a second metric g in which coincides with g in > R} but 

which is globally a small long range perturbation of the Minkowski metric. Let ip be the 
forward solution to Dgip = f. Then we set 



The estimate (1.5 I holds for the metric g, therefore we obtain 

sup E[^]{t) 



< 



Then the same bound holds as well for 0. Furthermore, we can compute the error 



f^{x>R-i)f + Pg,x>B]i' 



This has compact spatial support, and can be easily estimated in as in part (a) 



□ 



3.3. Proof of Theorem [3T2| , Given / G LE* we split it into 

/ Xehf Xpsf ~t" Xoof- 

For the last two terms we use part (b) of the Propositions 3.3|3.10 to produce approximate 



near the photon sphere, respectively near infinity. Adding them up 
we obtain an approximate solution 



solutions (t>ps and 



f. Due to ( 3.24| and (3.451 we obtain for 0o the bound 



sup£;[0o](*) + 110 



'Oil LB ^ 



l/ll 



for the equation DgCf 
(3.47) 

V 

while the error 

/l = ng(0ps + 0oo) - / 

is supported away from r = 3A/ and r = oo and satisfies 



I/: 



1\\LE, 



II/iIIl^< 11/11 



LE* 



Then we find = 0c, 

0i[O] = 0[O]-0o[O], 

By Theorem 1.2 we obtain the LEq bound for 0i. It remains to improve this to an LE 



01 by solving 

b.^heLE*, 



bound for 0i. By part (a) of Proposition 3.10 we can estimate ||xoo0||le 



Near the photon sphere we would like to apply part (a) of Proposition 3.3 to Xps*/*- 



However we cannot proceed in an identical manner because part (a) of Proposition 3.3 
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does not involve the Cauchy data of at i = 0, and instead applies to functions 4> defined 
on the full real axis in t. To address this issue we extend 0i backward in t to the set 
M'n, by solving the homogeneous problem Dg^i = in A^^, with matching Cauchy data 
on the common boundary of A^i? and Al'^. The extended function belongs to both 



LE{M.b) and LE{M'j^), and now we can estimate Xps^l^i via part (a) of Proposition 3.3 

4. StRICHARTZ ESTIMATES 

In this section we prove Theorem |1.4[ The theorem follows from the following two 
propositions. The first gives the result for the right hand side, /, in the dual local energy 
space: 

Proposition 4.11. Let (p,p,q) be a nonsharp Strichartz pair. Then for each £ LE 
with Ug4) e LE* + LiL^ 

we have 



(4.48) 



\LE-+LlL^ 



The second one allows us to use L^^L''^ in the right hand side of the wave equation. 

Proposition 4.12. There is a parametrix K for Dg so that for all nonsharp Strichartz 
pairs {pi,Pi,qi) and {p2,P2,q2) we have 

(4.49) supE[Kf]{v) + E[Kf]{E+) + + ||VA7|li.i^-p,.,, < ll/ll 



and the error estimate 
(4.50) 



\OgKf- fW^E'+LlL^ < ll/ll 



We first show how to use the propositions in order to prove the Theorem. 

Proof of Theorem \L4\ Suppose that Ugcf) — f with / e Lp'^Rp^'I'i, We write cj) as 

= </)! + L<f 



with K as in Proposition 4.12 By (4.49 1 the Kf term satisfies all the required estimates; 
therefore it remains to consider (pi. Using also (4.501 we obtain 



'gVl\\LE'+L\V 



+ £;[^i](o)<£;[0](o) + ll/ll 



iP2H''2.92 ' 



Then Theorem |3.2| combined with Duhamel's formula yields 



hWlE + ||ng0i|lL.+iii2 +supi?[0i](^) < i?[0](o) + ll/ll 



^P2^P2.<!2 ■ 



Finally the L^^H ''^''^ bound for V(/)i follows by Proposition 
We continue with the proofs of the two propositions. 



4.11 



□ 



Proof of Proposition 4-11 By Duhamel's formula and Theorem 3.2 we can neglect the 
L^L^ part of Og(j>. Hence in the sequel we assume that Dg^ e LE* . 

We use cutoffs to split the space into three regions, namely near the event horizon, 
near the photon sphere and near infinity, 

= Xeh(l> + Xps(f> + Xoo4'- 
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Due to the definition of the LE and LE* norms we have 

+ Eixpsm) + WXps^Ie,,, + \\Og{xps(t>)\\lE;^ 

+ E[xo.m) + iixoo^iiL^ + ipgixo^miEi,- 

Proving this requires commuting Dg with the cutoffs. However this is straightforward 
since the LE and LE* norms are equivalent to the H^, respectively L^, norm in the 
support of Vxe/i, Vxps and VXoo- 

It remains to prove the L^H^P'"^ bound for each of the three terms in V0. We consider 
the three cases separately: 

I. The estimate near the event horizon. This is the easiest case. Given (j) 
supported in {r < llAf/4}, we partition it on the unit scale with respect to v, 

(f>^^x{v- j)(l>, 

iez 

where x is a suitable smooth compactly supported bump function. Commuting the cutoffs 
with Og one easily obtains the square summability relation 

Y.\\^(^~^^'l'\\H^ + \\°aix{i~jmh+E[x{v-j)m<U^^^^^ 

where the energy term on the left is nonzero only for finitely many j. Since each of the 
functions x{v ~ j)4' have compact support, they satisfy the Strichartz estimates due to 
the local theory; see |32], |39], H?]. The above square summability with respect to j 
guarantees that the local estimates can be added up. 

II. The estimate near the photon sphere. For (j> supported in {5M/2 < r < 5M} 
we need to show that 

We use again the Regge- Wheeler coordinates. Then the operator Dg is replaced by Lfiw 
The potential V can be neglected due to the straightforward bound 

m\\LE;^<MLE,^. 

Indeed, for (j) at spherical frequency A we have 

\\V<P\\le;^, < |ln(2 + A)|5||0|U. < A|ln(2 + A)r5|10|U. < UUe,^. 

We introduce the auxiliary function 
By the definition of the LEpg norm we have 

(4.51) \mH^<\mLE,.- 

We also claim that 

(4.52) II WV'IU^ < II0IIle,„ + II W0IU^. 

Since A~} is L'^ bounded, this is a consequence of the commutator bound 

[Ap} , Ljiw] ■ LEps — > L'^, 
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or equivalently 

(4.53) [A~,\LRw]Aps: ^ L''. 

It suffices to consider the first term in the symbol calculus, as the remainder belongs to 
OPSf g, mapping to for all S > 0. The symbol of the first term is 

r*, A) = {a;/(A),e' + r~\r - 2M)X^}aps{X) 
and a-priori we have q e Sl'^^. For a better estimate we compute the Poisson bracket 

* U -la^ t 1 ^^ 4€r*-2g^,.(r-3(r-2M)) 
q(tr ,A) = ap, (A)7y(y,lnA) _^ ;^-2^2 

where y — r*^ + A^^^^. The first two factors on the right are bounded. The third is 
bounded by A since dr*{r~^{r — 2M)) vanishes at r* = 0. In addition, q is supported in 
1^1 < A. Hence we obtain q e XS^_g g. Then the commutator bound (4.531 follows. 



Given (4.511 and (4.521, we argue as in the first case, namely we localize -0 to time 
intervals of unit length and then apply the local Strichartz estimates. By summing over 
these strips we obtain 

l!Vi^||iPH-p., < UWle,^ + II W<^IIl2 

for all sharp Strichartz pairs {p,p,q). 
To return to we invert A~^, 

(j) = Aps^ + (1 - ApsA-})(l). 
The second term is much more regular, 

||V(1 - Ap,A~})^\\L2H^-^ < UWle,^, S > 0; 
therefore it satisfies all the Strichartz estimates simply by Sobolev embeddings. 

For the main term Apsip we take advantage of the fact that we only seek to prove 
the nonsharp Strichartz estimates for 0. The nonsharp Strichartz estimates for ip are 
obtained from the sharp ones via Sobolev embeddings, 

3 3 

||VV'||/i--P2,92 < ||Vl/;||ff-pi,,l, - +P2=- + Pi, Pi < P2- 

To obtain the nonsharp estimates for (f> instead, we need a slightly stronger form of the 
above bound, namely 

Lemma 4.13. Assume that I < qi < q2 < oo. Then 

3 3 

(4.54) \\ApsU\\H-P2.i2 < \\u\\h-pi.<ii, - +P2 = - + Pi- 

12 Ql 

Proof. We need to prove that the operator 

B = Op^ie + + l)-'tAp,Op'"{e + X^ + l)'^ 
maps L"^^ into L*^. The principal symbol of B is 

bo{r*,tX) = ie + A^ + l)'^aps{r*,^,X), 
and by the pdo calculus the remainder is easy to estimate. 
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The conclusion of the lemma will follow from the Hardy-Littlewood-Sobolev inequality if 
we prove a suitable pointwise bound on the kernel K of 6™, namely 

(4.55) |ifK,^i,r;,W2)| < (ki*-r*||c^i-c^2n"'+^"*. 

For fixed r* we consider a smooth dyadic partition of unity in frequency as follows: 

i = X{iei>A}+ x{^~t^}[xm<^o}+Yxm~v}\, 

/J dyadic \ i'=i'o / 

where i/q = I'd (A, r*) is given by 

In i^o(A, = In A + maxjln r*, — \/ln A}. 
This leads to a similar decomposition for 6oi namely 

&0 = &00 + X! I ^M,<i/o + X! ^^"^ I ■ 

In the region |f | > A the symbol 6o is of class S^^^^^ ^ which yields a kernel bound for 6oo 
of the form 

|i^oo(r-I,c.i,r;,a;2)| < {Vl - + 1^1-6.21)-^-"^+"^ 
The symbols of h^y are supported in {|^| m v, \ k, /i}, are smooth on the same scale and 
have size \\i{v~^ . Hence their kernels satisfy bounds of the form 

and similarly for Kfj,^^„g. Then ( |4.55| follows after summation. □ 
III. The estimate near infinity. 

Let us first recall the setup from [29]. We fix a Littlewood-Paley dyadic decomposition 
of frequency space in M.^, 

oo 

1= J2 suppsfeC{2'=-i < 1^1 <2'=+i}. 

A; — — OO 

Functions w in M x M"^ which are localized to frequency 2^ are measured in 

(4.56) \\u\\x, - 2'=/2||^.|U.(^^_^) + sup |||xri/2^|U.(^^), 

j>~k 

where 

= M X {2^ < \x\ < 2^+^}, A<j ^Rx {\x\ < 2^}. 
As in [29^ , by X° we denote the space of functions in M x K'^ with norm 



(4.57) Wuf^o^ W^'^^Wx. 

and by the dual norm 



k 

k — — OQ 



where Xj^ is the dual norm of X^^. 

One can establish the following (see |29l Lemma 1]) 
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Lemma 4.14. The following inequalities hold: 

(4.58) sup2-J'/2||Vu|U2(^^) <||Vm|Uo 

3 

and its dual 

(4.59) M\yo<\\u\\lei,. 

For small deviations from the Minkowski metric, one can also establish stronger local 
energy estimates involving the X° and F° norms; more precisely, one can prove (see |29l 
Theorem 4]): 

Lemma 4.15. Let g be a sufficiently small, long range perturbation of the Minkowski 
metric. Then, for all solutions u to the inhomogeneous problem DgU = / one has 

\\vu\\i^^. + iiv«ii^o < E[um + mi-.+L^^.. 

We now return to proving our estimate. For supported in {r > AM} we need to 
show that 

(4.60) \\^K.H-.. ^ ^['^KO) + UWIe., + \PMIei,- 

For large R we split (j) into a near and a far part 

^ = X>i?.0 + X<R(t> 



and estimate 



mm + uwIe,, + \PMiEi, > E[x>R4>m + \\x>rWle,. + \\^a{x>BMiEi, 

+ E[x<R^m + ||x<if^ll?fi + l|ng(x<i?0)lli- 

The term x<r4> has compact support in r and can be treated as in the first case (i.e. near 
the event horizon). Hence without any restriction in generality we can restrict ourselves 
to the case when (f) is supported in {r > R}. But in this region the operator Ug is a small 
long range perturbation of □; therefore the results of [29] apply. More precisely, from 
[291 Theorem 7(a)] we obtain 

I|V<^I1L^-,., < E[m + llV^lllo + |ln,0|lL.. 



This does not directly imply ( 4.60[l, s ince the X° norm is stronger than LEm- However, 



we can apply Lemma 4.15 and ( |4.59 l to obtain the bound 



\'^n\o<E[m) + \PMiEx 



□ 



Proof of Proposition \4^T2[ We split / into 

/ Xehf ^" Xpsf ^" Xc<if 

and construct the parametrix separately in the three regions. 

I. The parametrix near the event horizon. We further partition the term Xehf 
into unit intervals 

Xehf = ^X{V - j)Xehf 
j 
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with X supported in [—1,1], so that each component has compact support in the region 

-Dj = {ro < r < 11A///4, j - 2 <i < j + 2}. 
Let ipj be the forward solution to 

°sV'j = Xi.V - j)Xehf- 

Due to the local Strichartz estimates for variable coefficient wave equations, we obtain 
the uniform bounds 

W'^i^jWLPiH-Pi'OiiD,) + l|VV'il|L~L2(n,) + Wi^jU^L^D,) < \\x{v ^ j)Xehf\\ ^p'^ jjp2.<,'2 ■ 

Next we truncate 4'j using a cutoff function x{v — j, r) which is supported in Dj and 
equals 1 in the support of x{v—j)Xeh- Then the bound above also holds for the truncated 
functions — x{v — j, r)^jj, 

(4.61) ||V0j||iPiH-pi,,i + ||V(/)j||l~l2 + \\<l>3\\L^mD,) ^ \\X{V - j)Xehf\\^p'^^P2.i'2- 

In addition, 

°9<?i'j - Xiv- j)Xehf = - 

therefore 

(4.62) \\^g<l^o - X{i - i)Xehf\\L- < \\x{v-j)Xehf\\^p'2HP2.,'2- 

Finally, by energy estimates we also obtain a bound for the energy of 0j on the future 
space-like boundary of Dj at r = rg, 

(4.63) l|V0j||L2(£,^.n{,.=ro}) + ll0ilU2(D^.n{r=ro}) ^ \\x{v - j)Xehf \\ ^p'2 h»2.i'2 ■ 

To conclude we set 

Kehf ^^<t>j- 

Summing up the bounds (4.611, (4.621 and ( |4.63[ l for (jij we obtain the desired bounds 
for Keh^ namely 

respectively the error estimate 

W^gKehf - XehfWh^ < \\Xehf\\^p'^jjP2,,'2- 

II. The parametrix near the photon sphere. We work in the Regge- Wheeler 
coordinates. Arguing as in the previous case we produce a parametrix Kps with the 
property that, for each / supported in {5Af/2 + e < r < 5M — e}, the function Kpsf is 
supported in {5M/2 < r < 5M} and satisfies the bounds 

and the error estimate 

\\LRwKpsf - f\\L-2 < \\f\\^p'2HP2.l'2- 

Then we define the localized parametrix near the photon sphere Kp^ as 

E^psf — ^ps E^psXps^psiXpsf") 

with Xps = 1 in the support of Xps and slightly larger support. Then we show that Kps 
satisfies the required bounds. 
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We recall that (/C2-P2i 92) is a nonsharp Strichartz pair. Then by (4.54l we can write 

for some other Strichartz pair (^3,^3,(73) with = p2 and < q2. Since A^^ is 
bounded, from the above bounds for Kps we obtain 



By using (4.541 with Aps replaced by the weaker operator A^^ we also obtain the 
iPiij-Pi,9i bound for KpJ: 

IIVXps/lllpj^-Pi.,! < \\'^KpsXpsAps{Xpsf)\\'LPH-P'i ^ IIXps/ll^p^^P2.gJ, 1 

where {p,p, q) is another Strichartz pair with p — pi and q < qi. 

It remains to consider the error estimate, 
(4-64) \\LRwKpsf - XpsfWhE'+LlL-^ ^ IIXps/lliP^^P2-5^ ' 

for which we compute 

LBwKpsf — Xpsf — [Lnw , Ap}]KpsXpsAps{xpsf) 

+ Ap^{LRwKps - I)XpsAps{XpsI) 

+ {Ap}xpsAps - Xps){Xpsf)- 

We consider each term in the above decomposition. For the first term, due to the 
bound for K, we need the commutator bound 

[Lbw, Ap}] : — > LE* 

or equivalently 

Aps[LRw,Aps ] : H L , 
which is almost identical to (4.531 and is proved in the same manner. 

The bound for the second term is a direct consequence of the error bound for K. 

Finally, for the last term we know that {Ap]-Aps — I) £ OPS'^q^*; therefore using 
Sobolev embeddings we estimate 



i\s XpsAps - Xps)iXpsf)\\^p'2^^ < WXpsfW^p: 



2/^P2.<l2 ■ 



This concludes the proof of (4.64 1 since 



LP'^H^ C i^-ff 5 + L^H^ C LE* + L^L^. 

p^ 

III. The parametrix near infinity. 

We now consider the last component of /, namely Xoof- For some large R we separate 
it into two parts, 

Xoof = iXoo - X>R)f + X>Rf- 

The first part has compact support in r; therefore we can handle it as in the first case (i.e. 
near the event horizon), producing a parametrix K^. For the second part we modify 
the metric g for r < i? to a metric g which is a small, long-range perturbation of □. We 
let ipoo be the forward solution to 

DgV'oo = X>Rf- 
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We consider a second cutoff function x>R which is supported in r > R and equals 1 in 
the support of x>R- Then we define 

It remains to show that satisfies the appropriate bounds, 

snpE[K>^f]it) + \\K>^f\\lE,, + \\yK>"'f\\%,H-^,.,, < llx>i?/ll^.i^p,.,^ , 
respectively the error estimate 

W^gK^^ f - X>Rf\\LElj ^ \\X>Rf\\i^p'2HP2.i!i- 

These are easily obtained by applying the following lemma to ^poo- 

Lemma 4.16. Let f £ LP2 7jP2,92. Then the forward solution ip to Ogip = / satisfies the 
bound 



(4.65) 



It remains to prove the lemma. This largely follows from |29l Theorem 6], but there 
is an interesting technical issue that needs clarification. Precisely, |29l Theorem 6] shows 
that we have the bound 



(4.66) 



By Lemma[4.14| we are left with proving that 



(4.67) 



sup2 2 ||?/'||l2(^^) < ||/||^p^^P2,,^- 



We note that this does not follows from Lemma |4.14| this is a forbidden endpoint of 
the Hardy inequality in |29l Lemma 1(b)]. 



However, the bound (4.671 can still be obtained, although in a roundabout way. Pre- 



cisely, from (4.661 we have 
(4.68) 



supEmt)<\\ff^,. 



LP2HP2.'l2 



for the forward in time evolution, and similarly for the backward in time problem. 

On the other hand, a straightforward modification of the classical Morawetz estimates 
(see e.g. |27]) for the wave equation shows that the solutions to the homogeneous wave 
equation Ogip = satisfy 



(4.69) 



sup2-3^||V'||i2(^^) <i?[^](0). 



Denote by ltysH{t, s) the forward fundamental solution for Dg and by H{t, s) its back- 
ward extension to a solution to the homogeneous equation, DgH{t, s) = 0. Combining 
the bounds ( 4.68[ l and (4.691 shows that 

2 



sup2-3J 

3 



H{t,s)f{s)ds 



< 11/11' ' • '• 
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Since -p'^ < 2, by the Christ-Kiselev lemma |TT], it follows that 

2 



j 



which is exactly (4.671. 



Hit,s)fis)ds 



< 



l/ll 



□ 



5. The critical NLW 



In this section we prove Theorem 1.5 We first consider (1.181 in the compact region 



A4c- We denote by ip the solution to the homogeneous equation 

and by Tf the solution to the inhomogeneous problem 

□,(r/) = /, TF|So = 0, ifT/|So = 0. 
Then we can rewrite the nonlinear equation ( 1.18| in the form 
(5.70) <j> = ij±T{(l)^). 

We define Sobolev spaces in Aic by restricting to A4c functions in the same Sobolev 
space which are compactly supported in a larger open set. By the local Strichartz esti- 
mates we have 



and 



At the same time we have the multiplicative estimate 



110' 



5|| < 

14 ^ 



Then for small initial data we can use the contraction principle to solve (5.70 1 and obtain 
a solution (f) S H'i''^{Mc)- In addition, still by local Strichartz estimates, the solution 
(j) will have finite energy on any space-like surface, in particular on the forward and 
backward space- like boundary oi Mc- Thus we obtain 

E[^]{i:-^)<E[<j>]{^o). 



It remains to solve (1.181 in AA^ (and its other three symmetrical copies). Using the 
(■D,r, cj) coordinates in Mr we define -0 and T as above, but with Cauchy data on S^. 



By the global Strichartz estimates in Theorem |l.4[ for (s,p) as in the theorem we have 
and 

\\Tf\\LPH-.p(MR) ^ ll/lliii^- 



In particular we can take p = 5 which corresponds to s 
we have 



By Sobolev embeddings 
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therefore 



< 



3 

'//TO- 



Hence we can solve (5.70 1 using the contraction principle and obtain a solution (j) £ Hw^^. 
This implies that (j)^ g L^L^, which yields all of the other Strichartz estimates, as well 
as the energy bound on the forward boundary of Mr. This concludes the proof of 
the theorem. 



References 

[1] S. Alinhac: On the Morawetz-Keet-Smith-Sogge inequality for the wave equation on a curved 

background. Publ. Res. Inst. Math. Sci. 42 (2006), 705-720. 
[2] P. Blue and A. SofFer: Semilinear wave equations on the Schwarzschild manifold I: Local decay 

estimates. Adv. Differential Equations 8 (2003), 595-614. 
[3] P. Blue and A. Soffer: The wave equation on the Schwarzschild metric II: Local decay for the 

spin-2 Regge-Wheeler equation. J. Math. Phys. 46 (2005), 9pp. 
[4] P. Blue and A. Soffer: Errata for "Global existence and scatttering for the nonlinear Schrddinger 

equation on Schwarzschild manifolds", "Semilinear wave equations on the Schwarzschild manifold 

I: Local decay estimates", and "The wave equation on the Schwarzschild metric II: Local decay for 

the spin 2 Regge Wheeler equation", preprint. 
[5] P. Blue and A. Soffer: Phase space analysis on some black hole manifolds. J. Funct. Anal. 256 

(2009), no. 1, 1-90 . 

[6] P. Blue and A. Soffer: Improved decay rates with small regularity loss for the wave equation about 

a Schwarzschild black /ioie, http://xxx.lanl.gov, math/0612168 
[7] P. Blue and J. Sterbenz: Uniform decay of local energy and the semi-linear wave equation on 

Schwarzschild space. Comm. Math. Phys. 268 (2006), 481-504. 
[8] P. Brenner: On Lp — Lpi estimates for the wave- equation. Math. Z. 145 (1975), 251-254. 
[9] N. Burq: Smoothing effect for Schrddinger boundary value problems. Duke Math. J. 123 (2004), 

403-427. 

[10] N. Burq, P. Gerard, and N. T^vetkov: Strichartz inequalities and the nonlinear Schrddinger equa- 
tion on compact manifolds. Amer. J. Math. 126 (2004), 569-605. 

[11] M. Christ and A. Kiselev: Maximal functions associated to filtrations, 3. Funct. Anal. 179(2001), 
409-425. 

[12] H. Christianson: Dispersive estimates for manifolds with one trapped orbit, Comm. Partial Differ- 
ential Equations 33 (2008), p. 1147-1174. 

[13] Y. Colin de Verdiere and B. Parisse: Equilibre Instable en Regime Semi-classique: I - Concentra- 
tion Microlocale. Commun. PDF. 19, (1994), p. 1535-1563 

[14] M. Dafermos and I. Rodnianski: The red-shift effect and radiation decay on black hole spacetimes, 
arXiv:gr-qc/0512119vl. 

[15] M. Dafermos and I. Rodnianski: A note on energy currents and decay for the wave equation on a 

Schwarzschild background, arXiv:0710.0171vl. 
[16] M. Dafermos, I. Rodnianski: A proof of Price's law for the collapse of a self-gravitating scalar 

field. Invent. Math. 162 (2005), 381- 457. 
[17] J. Genibre and G. Velo: Generalized Strichartz inequalities for the wave equation. 3. Funct. Anal. 

133 (1995), 50-68. 

[18] S. W. Hawking and G. F. R. Ellis: The large scale structure of space-time. Cambridge Monographs 
on Mathematical Physics, No. 1. London, New York: Cambridge University Press 1973. 

[19] L. Kapitanski: Some generalizations of the Strichartz- Brenner inequality. Leningrad Math. J. 1 
(1990), 693-726. 

[20] M. Keel, H. Smith, and C. D. Sogge: Almost global existence for some semilinear wave equations. 

3. Anal. Math. 87 (2002), 265-279. 
[21] M. Keel and T. Tao: Endpoint Strichartz estimates. Amer. J. Math. 120 (1998), 955-980. 
[22] C. E. Kenig, G. Ponce and L. Vega: On the Zakharov and Zakharov-Schulman systems. 3. Funct. 

Anal. 127 (1995), 204-234. 



STRICHARTZ ESTIMATES ON SCHWARZSCHILD BLACK HOLE BACKGROUNDS 



43 



[23] J. Kronthaler Decay Rates for Spherical Scalar Waves in the Schwarzschild Geometry. 
arXiv: 6709. 3703 

[24] B. Kay and R. Wald: Linear stability of Schwarzschild under perturbations which are nonvanishing 
on the bifurcation 2-sphere Classical Quantum Gravity4 (1987), 893-898 

[25] I. Laba and A. SofFer: Global existence and scattering for the nonlinear Schrddinger equation on 
Schwarzschild manifolds. Helv. Phys. Acta 72 (1999), 274-294. 

[26] H. Lindblad and C. D. Sogge: On existence and scattering with minimal regularity for semilinear 
wave equations. J. Funct. Anal. 130 (1995), 357-426. 

[27] J. Metcalfe and C. D. Sogge: Long-time existence of quasilinear wave equations exterior to star- 
shaped obstacles via energy methods. SIAM J. Math. Anal. 38 (2006), 391-420. 

[28] J. Metcalfe and C. D. Sogge: Global existence of null-form wave equations in exterior domains. 
Math. Z. 256 (2007), 521-549. 

[29] J. Metcalfe and D. Tataru: Global parametrices and dispersive estimates for variable coefficient 
wave equations, arXiv:0707.1191 

[30] J. Metcalfe and D. Tataru: Global parametrices and dispersive estimates for variable coefficient 
wave equations in exterior domains, arXiv:0806.3409vl. 

[31] C. W. Misner, K. S. Thome, and J. A. Wheeler: Gravitation. San Francisco, CA: W. H. Freeman 
and Co., 1973. 

[32] G. Mockenhaupt, A. Seeger, and C. D. Sogge: Local smoothing of Fourier integral operators and 

Carleson-Sjolin estimates. J. Amer. Math. Soc. 6 (1993), 65-130. 
[33] C. Morawetz: Time decay for the nonlinear Klein-Gordon equations. Proc. Roy. Soc. Ser. A. 306 

(1968) , 291-296. 

[34] S. Nonnenmacher and M. Zworski: Quantum decay rates in chaotic scattering, preprint (2007). 
[35] H. Pecher: Nonlinear small data scattering for the wave and Klein-Gordon equations. Math. Z. 
185 (1984), 261-270. 

[36] R. Price: Nonspherical perturbations of relativistic gravitational collapse. L Scalar and gravita- 
tional perturbations Phys. Rev. D (3) 5 (1972), 2419-2438 
[37] J. V. Ralston: Solutions of the wave equation with localized energy. Comm. Pure Appl. Math. 22 

(1969) , 807-823. 

[38] W. Schlag, A. SofFer, W. Staubach: Decay for the wave and Schroedinger evolutions on manifolds 

with conical ends, Part IL arXiv:0801.2001 
[39] H. F. Smith: A parametrix construction for wave equations with C^'^ coefficients. Ann. Inst. 

Fourier (Grenoble) 48 (1998), 797-835. 
[40] H. F. Smith and C. D. Sogge: Global Strichartz estimates for nontrapping perturbations of the 

Laplacian. Comm. Partial Differential Equations 25 (2000), 2171-2183. 
[41] J. Sterbenz: Angular regularity and Strichartz estimates for the wave equation. With an appendix 

by I. Rodnianski. Int. Math. Res. Not. 2005, 187-231. 
[42] W. A. Strauss: Dispersal of waves vanishing on the boundary of an exterior domain. Comm. Pure 

Appl. Math. 28 (1975), 265-278. 
[43] R. S. Strichartz: A priori estimates for the wave equation and some applications. J. Funct. Anal. 

5 (1970), 218-235. 

[44] R. S. Strichartz: Restrictions of Fourier transforms to quadratic surfaces and decay of solutions 
of wave equations. Duke Math. J. 44 (1977), 705-714. 

[45] D. Tataru: Strichartz estimates for operators with nonsmooth coefficients and the nonlinear wave 
equation. Amer. J. Math. 122 (2000), 349-376. 

[46] D. Tataru: Strichartz estimates for second order hyperbolic operators with nonsmooth coefficients 
IL Amer. J. Math. 123 (2001), 385-423. 

[47] D. Tataru: Strichartz estimates for second order hyperbolic operators with nonsmooth coefficients 
in. J. Amer. Math. Soc. 15 (2002), 419-442. 

[48] D. Tataru: Parametrices and dispersive estimates for Schroedinger operators with variable coeffi- 
cients, to appear, Amer. J. Math. 

[49] F. Twainy: The Time Decay of Solutions to the Scalar Wave Equation in Schwarzschild Back- 
ground Thesis, University of California San Diego 1989 

[50] R. M. Wald: Note on the stability of the Schwarzschild metric. J. Math. Phys. 20 (1979), 105671058 

[51] R. M. Wald: General relativity. University of Chicago Press, Chicago, IL, 1984. 



44 JEREMY MARZUOLA, JASON METCALFE, DANIEL TATARU, AND MIHAI TOHANEANU 

Department of Applied Physics and Applied Mathematics, Columbia University, New 
York, NY 10027 

Department of Mathematics, University of North Carolina, Chapel Hill, NC 27599-3250 
Department of Mathematics, University of California, Berkeley, CA 94720-3840 
Department of Mathematics, University of California, Berkeley, CA 94720-3840 



